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Abstract 

Adopting the complete approach of the perturbative QCD (pQCD) and updated parton 
distribution functions, we have estimated the hadronic production of P-wave excited P c -states 
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(B*j L=1 ). In the estimate, special care on the relation of the production amplitude to the deriva- 



tive of wave function at origin of the potential model is payed. For experimental references, main 
uncertainties are discussed, and the total cross sections and the distributions of the production 
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with reasonable cuts at the energies of Tevatron and LHC are computed and presented. The 

0\ : 

results show that P-wave production may contribute to the P c -meson production indirectly by a 
factor about 0.5 of the direct production, and with such a big cross section, it is worth further to 



study the possibility to observe the P-wave production itself experimentally. 
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I. INTRODUCTION 

B c meson has been observed experimentally , and the observations are in consistency 
with the theoretical prediction. In addition to the experimental progress, as theoretical 
studies of the meson B c are getting deeply and widely 
H 16, 17, Q, 

interesting properties of the meson B C} such as its lifetime, mass, decay 
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branching ratios etc, are explored more and more precisely, the B c mesons can be used for 
generating tagged B s mesons (B c meson has very big decay branching ratio to B s meson) and 
for studying two heavy flavors c and b simultaneously etc, so B c physics is attracting more 
and more attentions. Considering the usages for experimental feasibility studies, recently we 
also wrote a generator named BCVEGPY [3] for the hadronic production of B c (S- wave), 
which is a Fortran program package and is in PYTHIA format 20]. 

Since excited states of B c , such as the P-wave states S*( 1 P 1 ) and P*( 3 P/)(J = 1,2,3) 
(for shortening, hereafter we use the symbol B* JL=1 to denote these four P-wave states) 
etc, may directly or indirectly (cascade way) decay to the ground state with almost 100% 
possibility via electromagnetic or hadronic interactions, thus the production of the excited 
B c states can be additional sources of the B c meson production, i.e. the B c mesons are 
produced 'indirectly' but promptly via excited state decay. If the indirect B c hadronic 
production and the direct one cannot be discriminated experimentally, but one would like 
to understand the production of the meson B c well, one certainly needs to know the indirect 
production fraction precisely. If one can discriminate the indirect B c production and the 
direct one and, furthermore, the excited B c states can be measured via their decay products 
exclusively with experimental techniques, then the potential models for the (cb) system may 



have various tests 



12j, therefore the theoretical estimate about the excited state production 



in advance certainly is very useful references for the experimental measurements. No matter 
the excited B c states can be measured exclusively or not, the theoretical estimate of the 
production of the excited states is requested. As known, it has a long story to understand 
the hadronic production of J/if>, even now there are still some problems e.g. experimental 
data still indicate not to agree with the theoretical prediction on the polarization of the 
hadronic produced J/ip and the discrepancies of the J/if> production in the B factory as 
well. While B c meson is flavored explicitly, so not only its hadronic production but also 
its excited states' hadronic production are simpler in comparison with the production of 



2 



the hidden flavored (cc) states. To understand the hadronic P-wave production not only 
is needed for understanding (cfe)-quarkonium (B c , B*, B* JL=1 , ■ ■ ■) production, but also will 
help to clarify up the situation about the hadronic production for hidden heavy flavor (cc) 
mesons. As a step, to study the production of the P-wave B c excited states and to write 
the generator accordingly are certainly interesting. 
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However, of the excited B c state hadronic production, only few authors of Refs. 
have studied the production of the P-wave ones B* JL=l so far. In Ref. 2l|, the hadronic 
B* JL=1 production is calculated with the fragmentation approach, which is comparatively 
simple and can reach to the leading logarithm (LL) level of perturbative QCD (pQCD). 
From the experiences for the S-wave production, the approach can reach to a high accuracy 
only when the production is in the region where the transverse momentum (p t ) of B c is 
very high ( Pt > 25 GevQ) to compare with the pQCD a\ approach. In addition, the 
complete pQCD a A s approach has a great advantage from the experimental point of view, 
that it retains the information about the c and b quark (jets) associated with the meson 
in the production. To retain the information about the c and b quark jets is more relevant 
experimentally, therefore, the complete pQCD oq approach to study the P-wave production 
of the B c excited states, as the case for .S'-wave production, is more favored, although it 



22J, |23| adopted the a s 



is the lowest order (LO) PQCD calculation. The authors of Refs. 
approach and found that the total cross-section is much greater than that predicted by the 
fragmentation approach, and the results from fragmentation approach can be compared with 
those from the approach only in the region where the transverse momentum of P* JL=1 , 
p t , is so high as p t > 30 GeV (higher than S'-wave production) j^ . 

The amplitude for the P-wave state production involves the derivative of the wave func- 
tion, so the derivation for the P-wave production is more complicated than that of the 
S'-wave state production a little. While, of the existence a* approach calculations 
the 'derivative' for the P-wave production is done numerically. 

Since the calculations are quite complicated, and the authors of Ref. j^J did the necessary 
derivative in the amplitude numerically, so to re-calculate and to verify the production 
calculations in terms of the approach is needed. To meet the needs, we also pay special 
attention on establishing the dependence correctly of the amplitude for the production on 
the derivative of the wave function at original which is fully determined in potential models. 
In fact, it is also to establish the correct relation between the matrix element appearing 
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in nonrelativistic QCD (NRQCD) |24| and the derivative of the wave function at original 
precisely. Therefore, we treat the hard subprocess amplitude in a strict way, that includes to 
write the amplitude under Bethe-Salpeter (BS) equation formulism J^, and to compute the 
derivative and to simplify the amplitude of the hard subprocess amplitude analytically as 
possible as we can (e.g. to find independent basic fermion strings and to expand the terms 
corresponding to Feynman diagrams accordingly) etc. Finally, the obtained amplitude before 
it to be squared is quite 'condensed'. To guarantee the Tightness of the calculations, we has 



done the comparisons between the results of ours and those in Refs. |22L |23( with the same 
input parameters, and done the checks of the gauge invariance for the amplitude numerically 
up to the computer abilities. 

The paper is organized: to follow Introduction, in Sec. II we highlight the dominant 
mechanism of gluon-gluon fusion, and present the basic formula for the B* JL=1 hadronic 
production in the complete ot 4 s approach. In Sec. Ill, we describe the input parameters 
for numerical calculations and present the results properly. Finally, in Sec. IV we present 
discussions and a short summary. 



II. FORMULATION AND TECHNIQUE 



Based on the factorization theorem of perturbative QCD (pQCD), the relevant hard 
subprocess plays a key role for hadronic production. In hadron collisions at high energies, 
for the production of the meson B c and its excited states, the gluon-gluon fusion subprocess 
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All of the discussions and 



gg — > B c (B*j L=1 ) + b + c is dominant 
calculations in the paper are based on the gluon-gluon fusion subprocess. 

At the lowest order of pQCD (af order), there are 36 Feynman diagrams for the gluon- 
gluon fusion process gg —> B* JL=1 + b + c totally, which can be schematically expressed as 
FIG^ so accordingly there are 36 terms in the production amplitude. As in Ref.|l2|, here we 
would like to establish the exact relation of each term of the amplitude to the wave function 
(derivative of the wave function at origin for P-wave production) of potential model (als o the 
relation between the matrix element in non-relativistic QCD (NRQCD) framework 241 and 
the wave function in potential model framework), thus we start with the 'rule' in Ref. |25j] to 
write down the explicit expression for the term of the amplitude corresponding to the kth 
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FIG. 1: The schematic diagram for the hadronic production of B* JL=1 , where the dashed box 
stands for the hard interaction kernel. k% and k 2 are two momenta for the initial gluons, q b i and 
q C 2 are the momenta for the outgoing b and c, P is the momentum of B* JL=1 . i, j, a and b are 
the color indexes of b, c and the corresponding gluons. 

Feynman diagram: 



M: 



A: 



C?««(to)i / ^ ( r ^ • X [(2S+1)p - J ^ ] (g) • T 2fc ) .(g c2 ) , (1) 

where X [(2S+1)Pj ' Jz \q) is the Bethe-Salpeter (BS) wave function for the bound state of (cb) in 
(2S+i)pj^ wnere (706 j s ^ ne co l or factor, i, j are the outgoing b and c quarks' color indexes 
and a, b are the color indexes of the initial two gluons respectively. stands for the 
structure of the k-th Feynman diagram that is between u(qbi) and the BS wave function 
2S+1 p Jt j z ] w hi cri includes the relative string of Dirac 7 matrices and the corresponding 
scalar part of the propagators, stands for the similar structure between the BS wave 
function 2S+1)p ^' J d(g) and v{q C 2). q b i and q C 2 are the momenta for the outgoing b quark 
and c antiquark. In the non-relativistic approximation, the BS wave function y^ S+1PjJz \q) 
can be written: 



Zr^WM* £ Z^L.^q).^.q). 
5 ~ A , Am b m c 

(qb2 - m b ) ■ (5s,o5 Sz fll5 + $s,i$s z ,\i£ x ') • (? c i + m c ) ■ (1A; SS Z \JJ Z ), (2) 

where P = Pr* , M = run* — m r + mj, are the momentum, the mass of B* JT _ A meson 
respectively (i.e. for shortening the notation, we adopt P and M to denote the momentum, 
the mass of B* JL=1 meson respectively) and q is the relative momentum between the two 
constituent quarks. Since q b2 and q c \ are the momenta of the b and c quarks inside the 
B* JL=1 , so they relate to the total and relative momenta of the bound state: 

n n m b m c 

q b2 = «iP -q, q cl = a 2 P + q, a x = ■ , a 2 = : • (3) 

m c + m b m c + m b 



For convenience in later usages, let us introduce 

and q\\ = \qjt\, hence accordingly we have d A q = dq\\d 3 qj_. The polarization vectors 

(6 A -P) = (e A '-P) = 0, A, A' = (1,0,-1); 

*&(q) stands for the 'P-wave scalar wave function', and (1A; SS Z \JJ Z ) is the Clebsch- Gordon 
coefficient for L-S coupling. In the paper, we will use a to denote the contraction between 
the Dirac 7 matrix and a momentum or polarization vector a, i.e. we use a instead of <ji. 
The spin structure of the BS wave function y^ 2S+1)Pj,Jz \q) defined in Eq.(j2j) is of the lowest 
order (up to 0(q)) for the P-wave state production. 

The P-wave scalar function ^(q) in Eq.(J2J), as a factor of BS wave function, relates to 
the derivative of the wave function at orig in ip' (0) in coordinate representation under the 
so-called instantaneous approximation [26] by the integration: 

dq\\d 3 q ± A f d 3 q± „ Q q± U , A „ , _ „.,a <*„,/, 



1 J ik^ qa * {q){€ ■ q) = l S (# S0( -^ )(e • ^ = ie V(0) • (4) 

„2 

Here q± is the unit vector ~^T ; an d 4>{~ stands for the P-wave 'scalar wave function' in 
the sense of potential model (the two quarks, as components of the bound state, must be in 
a relative space-like distance, hence <p should be a function of the variable jh)- Substituting 
Eq.Q into Eq.©, we obtain 

,2 



1 



^X l(28+1)P ^(q) ^ Y,^-^ X -^)^SS Z \JJ Z ) 

* ' S z ,\,\ 



• {—L=(-P + M) ■ (63,063,^5 + 5 s ,i5 Sz ,xie x ' 
k zv M 



- ( 



M 



X -Amf ) m c / 1 

+ ati-P + M) ■ (5 s ,o5 Sz ,ol5 + 5s,i5s z ,\>t X ')q±] + 0(q 2 ± )} 



,2 



S z ,\,\' 

,2 



where 



+ 0(qD), (5) 



A x ss z = -A^i-P + M ) ■ (*s,o*A',o7fi + hi^s z ,ye x ') , 



B 



a*f{-P + M) ■ (5<j,o<5a',q75 + S s ,i5 Sz ^e y ) 



V(t E _ 

\4m b m, 

+a 1 {-P + M) ■ (<5 s>0 «Ja',o75 + h,xh,^')l ,t 



X' u 

Note that at the lowest order relativistic approximation the terms proportional to B S g 
do not contribute to the production of the S- wave B C (B*) meson at all, however, to the 
production of the P-wave states B* JL=1 , they do contribute, so we have to keep the terms 
with care. To reach to the lowest order of the relativistic approximation, according to Eq. ([TJ 
the next step is to do the expansion of the and T 2k about g M up to 0(q 2 ) for the P-wave 
production, i.e. 

T lk = T° lk + T» k ■ q, + 0(q 2 ) , T 2k = T° 2k + T% k ■ q, + 0(q 2 ) , (6) 

where T^, r^, and Y 2k does not depend on g M (or q^) at all. Substituting the above 
equations into Eq.flTJ and carrying out the integration over d 4 q = dq\\d 3 q± with the help of 
Eq.(jlJ), we obtain 



M^=C^\0) £ (l\;SS z \JJ z )e*. 

S Z ,X,X' 

u s (q bl ) (i* ■ ■ T° 2k + T° lk ■ A$ 8a ■ T» k + T\ k ■ ■ T° 2k ) v s ,(q c2 ). (7) 

For a specific P-wave state, summing over the explicit A and A' for the Clebsch- Gordon 
coefficients, Eq.(J7j) can be further simplified as 

= c 2S+lp ^>(o)cz% 

■flW (r& ■ pf +1Pj]Jz ■ n k + r? fc ■ e? s+1 w* ■ r£ fc + r? fc • f^w- ■ r° fc ) v(q c2 ), (8) 

where the overall factor q (2S+1) Pj anc j ^he functions E^ s+1Pj ^ Jz and F^ s+1Pj ^ Jz are computed 
separately for different P-wave states B* JL=1 . 

For 1 Pi state, J = 1 , J z = A and S = 0, we obtain 

(9) 

Here and hereafter an extra factor 1 / \fW c will be implicitly included due to the fact that 
the meson is in color-singlet state. 
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The coefficients for the 3 Pj(J = 0, 1, 2) states with S = 1 can be obtained with the help 
of the relations 0, : 

£<1A; lA'|00)eJe? = ^ - <? M „) , ( J = 0) , (10) 



A, A' 



£<lA;lA'|lJ,)e£ei' = iJl e ^^ef , (J = 1,J Z = -1, 0, 1) , (11) 



A, A' 

£<1A; 1A'|2 J 2 )e^' = , {J = 2,J Z = -2, -1, 0, 1, 2) , (12) 

A, A' 

and the polarization vector and tensor e^ v {J z ) obey the projection relations 

= (^-^)=*V, (is) 

^1 6 fiv^ pa = l^Plip^VQ + 'PupP^ ~ -^Pta/PpQ i (14) 

for 3 P state (A = 0): 



E [ 3 Po]A = ( _p + M)(M +p) F [3Po]A = _ 3[(^-a 1 )P + M] 3 p Q = 1 

(15) 

for 3 P X state (A = -1,0,1): 
M 3 ^ = ,(P + M)e/^ 7p P^, F ?P^ = _ P ^-Me^ 3 Pl = 1 

v 'pquj aia2 , 2v^at;m 3 / 2 v ; 

Finally for 3 P 2 state (J, = -2, -1, • • • , 2): 

4 3P2,Jz = (^ + %74. P^=0, = (17) 

If one calculates cross sections of the production by summing all the 36 terms (cor- 
responding to the Feynman diagrams), and then having the result squared directly, the 
squared amplitude will be too long to deal with, even in terms of computer, and is really 
time-consuming to reach to the final results for the cross sections. Instead of summing the 
terms directly, we adopt the FDC package [2^ to generate the Fortran program for each 
term and develop a technique in FDC further. Namely we firstly establish a complete set 
of ' basic spinor lines', such as that they are constructed by the multiplication of Dirac 
7 matrixes in a certain way as "px^Vz ■ • ■, then we consider them as bases to expand every 
terms of the amplitude and sum up all the terms according to the expansion (the coefficients 
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of the 'bases' are summed respectively). In this way, the amplitude and its square become 
quite condense, thus the efficiency for computing the amplitude squared is raised greatly. 

In order to simplify the amplitude manipulation under Fortran program, it is better to 
adopting the linear polarizations instead of the circular ones because the linear ones are 
real numbers. Thus we do so, and for convenience we put the explicit form of the linear 
polarization vector e M and that of tensor e^ v in Appendix I, which are used in our calculations 
when summing up the final polarizations of B* L+1 . 

The color factors for the terms corresponding to the Feynman diagrams should be also 
treated correctly. For the gluon-gluon fusion subprocess, we choose the three independent 
color factors |37( as follows: 



3\/26 

= ^7= (2(T°T% + 2(T b T% - 35 tJ Tr[T a T b ]) , (18) 

TT = ^^(-ll(T a T% + 2{T b T a ) l3 - 35 l3 Tr[T a T b }) , (19) 

T% ab = -^L(-3(T b T% - 6 tJ Tr[T a T b }) , (20) 



with the normalization: (T™' a *T"' a ) = 5 mn (m = 1, 2, 3). All the color factors C^ - in 

a,b,i,j 

Eq.(fTf can be expanded as the linear combination of these three independent color factors. 
Therefore the amplitude can be grouped into three terms by 'color', 

m s,jj z = Mf' JJz = {T™' ab M% JJ *). (21) 

fc=l,---,36 m=l,2,3 

For various un-polarized cross sections of the subprocess (in given J and S), we have the 
square of the amplitude: 

\M S > J \ 2 = |M 5 ' JJz | 2 . (22) 

According to pQCD factorization theorem, the hadronic production cross section is for- 
mulated as 

da = Y^[d x i[ dx 2 F l Hi (x u /j, 2 f ) x F J H2 (x 2 , n 2 F )da ij ^Bi JL ^ 1 x{.xi, x 2 , n%, /i 2 , Q 2 ) , (23) 

u 

where F l H ( x, fJ?jr) j Fjj 2 (x , /if?) are the parton distribution functions (PDFs) of the i and j 
partons in the hadrons Hi, H 2 respectively, fi 2 is the 'energy scale squared' where renormal- 
ization for the subprocess is made; Q 2 is the 'characteristic energy scale of the subprocess 
squared' i.e. when setting fi 2 = Q 2 , then LO pQCD calculations for the subprocess can 



obtain the best results, thus without emphasis fi 2 = Q 2 is always set; and fi 2 F is the 'energy 
scale squared' where the factorization about the PDFs and the hard subprocess is made. 
Usually for LO to obtain the best results, the factorization and 'renormalization' are carried 
out at the same energy scale i.e. \j? F = Q 2 , thus later on we take /i 2 = fip = Q 2 and define 
d&ij-+B* j ^xi^i, X2, Q 2 ) = daij^B*j L _ 1 x(^i, X2, [i F , /i 2 , Q 2 ) except one case when estimating 
the uncertainty from LO and the ambiguity of the choices about [i 2 F and /i 2 = Q 2 . 

Having all the above preparations, we base on Eq. (}2*3*|) technically to write a program for 
numerically computing the hadronic production cross sections at various energies (31 1. 



III. NUMERICAL RESULTS 



Before doing numerical calculations, we have done the necessary checks. First of all, 
we numerically check the gauge invariance. We find the gauge invariance is guaranteed 
at the used computer ability level for our formula and computer program to compute the 
amplitude. Then we further compute the cross sections numerically for the subprocess 



g + g — > B*j L=l + b + c with the input parameters as taken in Ref. 



numerical results as those in Ref. 



23, and the same 



2J are obtained (see TABLE [fl FigOalso Ref. [22]). Thus, 



we confirm the results of Ref.|22|] for the subprocess, and are quite sure the correctness of 
our formula and program. 

We should note here that the approximation Eq.(J2J), which is also taken in Ref. and 
is similar in the cases for (cc) and {bb), the heavy flavor hidden systems, indicates that 
the equations P = q c \ + qm-, <?ci = m c' Qb2 = m l an d P 2 — m B* must be satisfied 

cJ,L — l 

simultaneously, furthermore, only when the equations are satisfied, the gauge invariance of 
the amplitude is guaranteed. Therefore, under the approximation, there are two possibilities: 
one is to take the same quark mass values in the P-wave and 5-wave state production, 
but one have to ignore the mass difference between the P-wave and S"-wave states i.e. 
rriB c = m B* JL _ 1 ] the other one is to consider the mass differences between the P-wave and 
S'-wave states i.e. ^ m B* JL _ 1 i n the P-wave production, but one have to take different 
mass values of c and b quarks from those in S'-wave production accordingly. The authors 
of Ref. took the later i.e. they took nib = 5.0 GeV, m c = 1.7 GeV for the P-wave state 
production in order to have a mass of B* JL=l about 6.7 GeV, although they took m# c = 6.3 
GeV, m c = 1.5 GeV and rrib = 4.8 GeV for S'-wave production in the meantime. As known 
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TABLE I: The total cross section for the hard subprocess g + g — > B* JL=1 + b + c (gluon fusion 
into P-wave excited states B* JL=l ) at different CM. energies. The input parameters are taken as 

n 

those used in Ref.[22j: m;, = 5.0Gel/, m c = 1.7GeV, M = 6.7GeV and the running a s is fixed to 
0.2 etc. 



CM. energy (GeV) 


20GeV 


40GeV 


60GeV 


80GeV 


lOOGeV 


200GeV 


a( 1 P 1 )(pb) 


0.184 


0.743 


0.657 


0.538 


0.439 


0.195 


a( 3 P )(pb) 


0.367 


0.207 


0.175 


0.141 


0.114 


0.0496 


a( i P 1 )(pb) 


0.346 


0.598 


0.503 


0.402 


0.324 


0.139 


a( 3 P 2 )(pb) 


0.721 


1.49 


1.31 


1.06 


0.862 


0.374 



that the cross section of the jet production g + g — > c + b + b + c is quite sensitive to the 
chosen values of the masses mb, m c , thus, we think that the jet production should not be 
affected by the differences for the P-wave and S'-wave production, and in fact that is also 
the requirement for the NRQCD formulism, of the two possibilities which cannot be avoided 
under the approximation Eq.(J2J), the former, i.e. to take w^* Ji _ 1 = % an d to keep the 
values of b and c quark masses to be the same as those for the S'-wave production, should 
be more relevant comparatively. Hence, in the paper we carry out all numerical calculations 
always to take the 'choice': the mass values of c, 6-quarks as m c = 1.5 GeV = 4.9 GeV, 
i.e., the same as those taken in S-wave production, so rrtB* JL _ 1 = vtlb c = 6.4 GeV, except 
the case when showing the uncertainties from the two possible choices. 

n 

As pointed out in Ref. 32J, there are quit a lot of uncertainty sources in the hadronic 
production for the S- wave B C (B*) meson, such as the variations about a s -running, the vari- 
ations in choices of the factorization energy scale, the adopted PDFs in different versions 
and various input parameter values relating the bound state etc.. For instance for definite- 
ness in check of our formula and computer program by comparison with the results of [22], 
we take pQCD coupling constant a s = 0.2, the derivative of the wave function at origin 
|P'(0)| 2 = 0.201GeV 5 , and the mass values of c-quark, 6-quarks and B* JL=1 are fixed as 
m b = 5.0GeV, m c = 1.7GeV, M = 6.7GeV etc i.e. all are fixed as those in for the 
subprocess gg — > B* JL=1 + b + c. 

To compute the production, we need to know the PDFs. Several groups, such as CTEQ 
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FIG. 2: B* L=1 — p t and B* L=1 — y differential distributions for the subprocess. The dashed line 
(up- middle) , the diamond line (bottom), the dotted line (down-middle) and the solid line (top) are 
for P%, 3 Po) 3 Pi an d 3 P2, respectively. The parameters are taken to be the same ones that have 
been used in Ref.[22|. 

Q, GRV Q and MRS Q etc., devote themselves to offer accurate PDFs to the world. 
By taking different PDFs, one may find that the results caused by different PDFs can be 
about 10%. Later on for definiteness, we will take CTEQ6L[33] and the leading order a s 
running with = 0.326 GeV. 



As pointed out in Ref. [32], at LO pQCD of all the uncertainties, that caused by the 
choice of the energy scale Q 2 for the renormalization of the subprocess (or say that for the 
factorization) is the largest in the hadronic production estimate, and for 5-wave production 
it may cause such uncertainty as great around a factor 30%. In fact, for P-wave production 
it is also the case. To see the uncertainty caused by Q 2 choice for the P-wave production, 
we try to calculate the P-wave production under the condition that the other uncertainty 
factors are fixed. As in Ref. 3^|, we take four types of Q 2 choices and put the results for 
total cross section in TablTTl respectively. The four types are: 

Type A: Q 2 = s , the squared CM. energy of the subprocess; 

Type B: Q 2 = M 2 = p 2 + m 2 B * , the squared transverse mass of the B* JL=1 meson; 

cJ,L — 1 ~ ' 

Type C, Q 2 = m 2 t = p 2 b + m 2 , the squared transverse mass of the b quark. 
Type D: Q 2 = Am 2 . 

From TablTTl about the total cross sections for the hadronic production of B* JL=1 at 
LHC and TEVATRON with the different types of the choices, one may observe that the 
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TABLE II: The total cross section (in unit nb) for the hadronic production of P-wave B c meson 
(B* JL=1 ) at LHC and TEVATRON with different types of factorization energy scales (A, B, C, D). 
Here m?, = 4.9 GeV, m c = 1.5 GeV and m^* = m c + rrih. 

cJ,L=l 





LHC (VS = 14. TeV) 


TEVATRON (VS = 1.96 TeV) 


Q 2 


A 


B 


C 


D 


A 


B 


C 


D 


<r{ X Pi){nb) 


4.738 


9.123 


9.825 


8.379 


0.2555 


0.6545 


0.7547 


0.5507 


a( 3 P )(nb) 


1.910 


3.288 


3.523 


3.036 


0.1161 


0.2563 


0.2966 


0.2149 


a{ 3 Pi){nb) 


4.117 


7.382 


7.304 


6.682 


0.2289 


0.5597 


0.6490 


0.4780 


a( 3 P 2 )(nb) 


10.18 


20.40 


21.71 


18.26 


0.5096 


1.350 


1.515 


1.102 



Q 2 dependence in the P-wave B* JL=1 production is much more stronger than that in the 
S-wave production. With the different choices, the total cross sections of the production 
can be varied so big as a factor 2.0 ~ 3.0. 

For definiteness in the rest estimates on the production we will adopt just Type B for 
Q 2 choice. 

For experimental references, we also calculate the distributions of p t and y (the B* JL=1 
transverse momentum and rapidity) precisely and draw the curves on the distributions in 
FigEJfor LHC and in Figgjfor TEVATRON. 

For comparison, in Figs l3l4l we also show the p t and y differential distributions for the 
production of the S-wave states 1 S'o and 3 Si. One may observe that the p t and y dependence 
of the differential distributions behaves quite similar to those of the 1 Sq state production. 
The summed up P-wave production cross section of the P-waves 1 Pi, 3 Po, 3 Pi, 3 P2 is smaller 
than the S-wave production, but it can reach to a fraction about 50% of the 1 S production 
cross section. 

From FIGs loTH one may also see that the p t and y differential distributions for the 
hadronic production of B* JL=1 are quite different at LHC and at TEVATRON due to the 
different CM. energies of the two colliders. 

To see the uncertainties from Q 2 choice, instead of variation on the choices with Q 2 = 
fip, the authors in literature, such as Ref.j^, also try Q 2 7^ p 2 F in Eq.(|23|) and see the 
uncertainty. Here following them, we calculate the summed pt- and ^/-distributions of all 
the P-wave P Ci l=i production for LHC and TEVATRON and present the uncertainty on 
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FIG. 3: The differential distributions of pt and y for the hadronic production of B* JL=1 at LHC. 
The dotted line (down- middle) , the diamond line (bottom), the dashed line (next to bottom), the 
solid line (up-middle) are for Pi, s Po, 3 Pi, 3 P2, respectively. For comparison, the results for 1 Sq 
and 3 S*i wave states are shown in triangle line (next to top) and dash-dot line (top) correspondingly 



(with parameters: m& = 4.9 GeV, m c = 1.5 GeV and niB c = ms* 



m c + m b ). 



the distributions in FIGs l5l6l respectively. One may see the uncertainty can be so great as a 
factor 4.0. 

In high energy hadronic colliders, the events with a small p t and/or a large rapidity y 
cannot be detected by detectors directly, so for experimental studies and for precise purposes 
in the estimates, only the events with proper kinematic cuts on p t and y are taken into 
account. In Ref. 2, we did the studies on the 5-wave production, thus for experimental 
references here we also try various cuts for the P-wave production. 

To study the cut effects on the production, here we consider the production in summing 
up all of the P-wave ( 1 Po, 3 Po, 3 Pi, 3 P2) production. Considering the abilities in measuring 
rapidity of B c for the detectors CDF, DO and BTeV at TEVATRON, and ATLAS, CMS 
and LHC-B at LHC, we try the rapidity cuts y cut ~ 1.5 and higher for p t distribution and 
the ptcut ~ 5 GeV and higher for y distribution. The results with the four cuts: y cut = 
(0.5,1.0,1.5,2.0) are put in FIG0 and the results with the four cuts: p tcut = 5,20,35,50 
GeV are put in FIGJBJ 

From FigQ we can see that the dependence of the p t distributions on y tcut for LHC 
is stronger than that for TEVATRON. This is because that at TEVATRON, the sizable 
rapidity distributions cover a smaller region in y, so generally speaking the rapidity cuts 
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FIG. 4: The differential distributions of pt and y for the hadronic production of B* JL=1 at TEVA- 
TRON. The dotted line (down-middle) , the diamond line (bottom), the dashed line (next to bot- 
tom), the solid line (up-middle) are for Pi, 3 Po, 3 -Pi, 3 P2, respectively. For comparison, the 
results for 1 Sq and 3 Si wave states are shown in triangle line (next to top) and dash-dot line (top) 
correspondingly (with parameters: mj = 4.9 GeV, m c = 1.5 GeV and tob c = ^tib* 3 h _ x = m c + mb). 

affect the production stronger than that at LHC. The correlations between p t and y can be 
seen more clearly from FiglHJ the y-distribution various pt-cuts. One may observe that the 
dependence of the differential distribution on rapidity y with different p tcu t is to exhibit a 
broader profile at LHC than at TEVATRON. 

We also would like specially to show the consequences due to the approximation Eq.Q 
with one by the two possibilities described above, i.e. the one taken by ourselves: mg« _ = 

B c = 6.4 GeV, m c = 1.5 GeV and mj = 4.9 GeV, and the other one taken by Refs.|3, 
mB h L-i = GeV, m c = 1.7 GeV and nib = 5.0 GeV quantitatively, because it is 
'fresh' in P-wave production. We calculate p t distributions for the P-wave ( 1 Pi, 3 Po, 3 Pi, 3 P2) 
production with the two possibilities respectively and draw the curves of the differential cross 
sections versus p t in FIGE1 It is obvious that with our present parameters for the constitute 
quarks, the distributions (hence the total cross-sections) are about two times bigger than 
that of the second choice for the constitute quark masses. Namely, the results mean that 
the consequences due to the approximation Eq.(j2J) are sizable and indicate improvement on 
the approximation is requested. 
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FIG. 5: The summed pt~ and y-distributions of all the P-wave P Cj l=i for different values of 
the factorization scale \j? f and the renormalization scale Q 2 at LHC. The upper edge of the 
band corresponds to {jl 2 f = AM 2 , Q 2 = M 2 /A and the lower edge corresponds to that of 
H F = M 2 /A, Q 2 = AM 2 . The solid line, the dotted line and the dashed line corresponds to 
that of n 2 F = Q 2 = M t , fJ 2 F = Q 2 = AM t , fj? F = Q 2 = M?/A respectively. 




dy p,(6eV) 



FIG. 6: The summed up pt- and ^-distributions of all the P-wave B c ^l=i for different values of 
the factorization scale up and the renormalization scale Q 2 at TEVATRON. The upper edge of 
the band corresponds to \j? f = AM 2 , Q 2 = M 2 /A and the lower edge corresponds to that of 
H F = M 2 /A, Q 2 = AM 2 . The solid line, the dotted line and the dashed line corresponds to that of 
H 2 F = Q 2 = M t 2 , fj 2 F = Q 2 = AMI, Vf = Q 2 = M t/ 4 respectively. 
IV. DISCUSSIONS AND SUMMARY 

In the paper, we have studied the hadronic production of the P-wave B* JL=l in a way as 
analytical as possible, i.e., the amplitude of the hard subprocess g + g — > B* L=1 + b + c is 
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FIG. 7: The summed up differential distributions for all the considered P-wave state production 
versus pt of B* JL=1 with a y cut (y C ut) at LHC (left) and at TEVATRON (right). Dashed line 
(next to top) with y cut = 2.0; dash-dot line (middle) with y cut = 1.5; dotted line (next bottom) 
with y cut = 1.0; diamond line (bottom) with y cut = 0.5 and solid line (top) without y cu t- 
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e 

T3 



T3 



-3 -2 -1 




FIG. 8: The summed up differential distributions for all the considered -P-wave state produc- 
tion versus rapidity y of B* JL=1 with a transverse momentum cut (ptcut) at LHC (left) and at 
TEVATRON (right). Diamond line (next to top) with ptcut = 5 GeV; dash-dot line (middle) with 
Ptcut = 20 GeV; dotted line (next to bottom) with ptcut = 35 GeV; dashed line (bottom) with 
Ptcut = 50 GeV and solid line (top) with ptcut = 0. 

reduced with the help of the FDC program and the technique to expand each term of the 
amplitude by the independent 'bases' for the 'elementary fermion strings' analytically then 
to sum up all of the 36 terms (corresponding to the 36 Feynman diagrams). Thus having the 
amplitude in expansion of the bases squared, we have made the numerical calculations very 
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FIG. 9: The p t distributions of all the P-wave B* L=1 ( 1 P 1 and 3 Pj (J = 1,2,3)) production 
at LHC (left) and TEVATRON (right). Of the same type of lines, the upper is for the choice: 
nib = 4.9 GeV, m c = 1.5 GeV, while the lower is for the choice: mj = 5.0 GeV, m c = 1.7 GeV, 
respectively. For LHC, we take \y\ < 1.5 and for TEVATRON, we take \y\ < 0.6. 



efficiently. The correctness for the amplitude and the program is tested by demonstrating its 
gauge invariance numerically up to the computer abilities and comparing the results for the 
subprocess with those in Ref . |22j by taking the same parameters. Having done the tests, we 
will add the whole program into the B c meson generator as the new version BCVEGPY2.0 
soonloT]. 

From the figures FIGs 12"IqT41 and the tables Tabs JIIHl one may see that, of the hadronic 
P-wave production, that of 3 P2 is the biggest and that of 3 Pq is the smallest. From the 
table TablHl and the figures FIGsEEl one ma Y see that, at LO, the uncertainties from the 
factorization and renormalization energy scale choice, i.e. on Q 2 ,^ choices for the P-wave 
production, are quite great, even greater than that of S'-wave production. Furthermore, from 
the figure FIGJH1 one may see that the approximation Eq.(J2J) causes a great uncertainty in 
the production estimate. 



Relating to the approximation, our results with the choice rag. 



m Bc = 6.4 GeV, 



m. 



Ref. 



1.5 GeV and m;, = 4.9 GeV can be almost two times bigger than the results obtained in 



22| with the choice tub* 



6.7 GeV, m c = 1.7 GeV and m b = 5.0 GeV at LHC and 



TEVATRON. The fact about the factor 2 can be understandable: as known, smaller quark 
masses will cause the cross section of the four quark production gg — > c + c + b + b bigger 
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and they will also cause an increase effect in the binding c + b — > B* JL=1 , thus when we take 
smaller masses of quarks and also the bound state mass M = mh + w, c in the meantime, 
we obtain such bigger results by the factor 2 than those in Ref.|22j. We think that for the 
production the true values (without the approximation Eq. (0) for a more accurate estimate) 
should be smaller than those we present here, due to a comparatively small value of the 
masses of B* JL=1 taken here, but the true values should be bigger than those obtained in 



Ref. 



22| due to there comparatively big values for the masses of c and b quarks taken for 



the factor of the four quark production inside the P-wave production respectively Since 
the factor of the involved production gg — > bbcc is guaranteed in gauge invariance, here 



221 not 



only the factor cb — > B* JL=1 is not well treated in the present results, while in Ref. 
only the factor cb — ► B* JL=1 but also the factor of the involved four quark production are 
not well-treated by comparatively bigger values for the quarks and the bound states, so we 
suspect that the true values for the production might be closer to the presented ones here, 
rather than those in Ref. (2^. Note that the term B^' s ^q^ in Eq.(JHJ) does not play any role 
for 5-wave B c and B* production at all, thus for the .S-wave production the approximate 
choices about the values of P = q c \ + q^ , m c , nib and m Bc ^ B *^ are not so sensitive as the 
P-wave production. 

The summed total cross section for all of the considered hadronic P-wave production 
can be so big as a half of the direct production of the ground state B c ( 1 Sq). Considering 
the the fact that almost all of the P-excited states B*( 3 S\) decay to the ground state B c , 
the contribution from the B*( 3 Si) production to the 'final P c -wave production', we may 
conclude that the contribution can be about 20% in total. From the pt and y distributions, 
one may see that P-wave and S-w&ve production behave quite similar, that no matter what 
Ptcut and y cu t are chosen, the P-wave production itself with such big cross sections is worth 
while to study the possibility directly measuring the P-wave P c -states seriously. Especially, 
the properties of the P-wave mesons B^ c L=1 ^ are crucial in understanding the mass spectrum 
of the c6-quarkonium states and testing the potential models. 
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APPENDIX A: THE LINEAR POLARIZATION VECTOR AND TENSOR 

As stated in the text, the linear polarization vector (tensor) contains real number only if 
the meson (bound state) is a vector (tensor) meson, so the linear ones are better than com- 
plex circular polarization ones in numerical calculations, thus the linear explicit expressions 
for them are quite useful in writing numerical programs. Hence, we present the precise ex- 
pressions for the linear polarization vector if the meson is a vector one, and the polarization 
tensor if the meson is a tensor one in the appendix. 

We set a coordinate first, in which z-axial is in the direction of the incoming hadron, 
and the momentum of the concerned meson (bound state) in the direction: P M = 
(Po, I P I sin 9 cos 0, |P| sin # sin 0, |P|cos#) with |P| = J P 2 — M 2 , i.e., 9 is the polar angle 
and is the azimuth angle in the coordinate. The space-like linear polarization vector for 
the meson can be expressed as: 

e£(P) : (0, cos # cos 0, cos 9 sin 0, -sin#) (Al) 
e£(P) : (0, -sin0, cos0, 0) (A2) 
<£(P) : i(|-P|, Posin#cos0, P o sin0sin0, P cos9) . (A3) 
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As the request of the polarization vector definition, they do satisfy the conditions: 

€i ■ P = , €i ■ €j = -Sij , (i, j = x, y, z) . (A4) 

Accordingly, the 'linear polarization tensor' e° (P) for B*(^P-2) can be constructed in 
the following way i.e. the five components (J 2 = 1,2, •• - 5) of the polarization tensor are 

— * — * — * 

constructed in terms of e x (P), e y {P) and e 2 (P): 

ja> _ J_/ u v _|_ e ti e v\ J"' — J^ffUfV i 

fc i ~ ^l^^Ttye^j, e 2 — ^^tj.t^j -t- t 2 t,j , 

er = 75(^ + ^-2e^), (A5) 

where again /x and z/ are the Lorentz vector indexes. Again, one may easily check that the 
polarization tensor components satisfy the conditions: 

e\.€j = -Sn , = g, u e7 = , P M • ef = • P^ = for (i, j = 1, • • • , 5). (A6) 
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